Measurement of Beam-Spin Asymmetries for 
7T + Electroproduction Above the Baryon Resonance Region 
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We report the first evidence for a non-zero beam-spin azimuthal asymmetry in the electropro- 
duction of positive pions in the deep-inelastic kinematic region. Data for the reaction ep — > c'tv + X 
have been obtained using a polarized electron beam of 4.3 GeV with the CEBAF Large Acceptance 
Spectrometer at the Thomas Jefferson National Accelerator Facility (JLab). The amplitude of the 
sin</> modulation increases with the momentum of the pion relative to the virtual photon, z. In the 
range z =0.5-0.8 the average amplitude is 0.038 ± 0.005 ± 0.003 for a missing mass Mx >1.1 GeV 
and 0.037 ± 0.007 ± 0.004 for Mx >1.4 GeV. 

PACS numbers: 13.60.-r; 13.87.Fh; 13.88.+e; 14.20.Dh; 24.85.+p 



The origin of the spin of the proton has become a 
topic of considerable experimental and theoretical inter- 
est since the EMC pj measurements implied that quark 
helicities account for only a small fraction of the nucleon 
spin. As a consequence, the study of the gluon polar- 
ization and the orbital angular momentum of partons 
have become of central interest. Single-Spin Asymme- 
tries (SSAs), measured in hadronic reactions for decades 
H @| , have emerged as important observables to access 
transverse momentum distributions of partons and the 
orbital angular momentum of quarks in the nucleon. 

In this paper we present the first measurement of a 
non-zero beam-spin asymmetry in the electroproduction 
of positive pions in deep-inelastic kinematic region. Re- 
cently measurements of target SSAs have been reported 
for pion production in semi-inclusive deep-inelastic scat- 
tering (SIDIS) by the HERMES collaboration 0, H 
for longitudinally polarized targets, and by the SMC col- 
laboration for a transversely polarized target Q. Such 
SSAs require a correlation between the spin direction of 
a particle and the orientation of the production (or scat- 
tering) plane, and have been linked to the orbital an- 
gular momentum of partons in the nucleon |8|, |9| . Two 
fundamental QCD mechanisms giving rise to single spin 
asymmetries were identified. First the Collins mecha- 
nism |lOl ITU El ] , where the asymmetry is generated in 
the fragmentation of transversely polarized quarks, and 
second the Sivers mechanism El III El El El El El , 
where it arises due to final state interactions at the distri- 
bution function level. The interference of wavefunctions 
with different orbital angular momentum, which is re- 
quired to generate the SSA 0, 0, 0, 0] , also yields 



the helicity-flip Generalized Parton Distribution (GPD) 
E [a li| that enters Deeply Virtual Compton Scatter- 
ing (DVCS) jHIIil and the Pauli form factor F 2 . The 
connection of SSAs and GPDs has also been discussed in 
terms of the transverse distribution of quarks in nucleons 

mm 

Physical observables accessible in SSAs include novel 
distribution and fragmentation functions such as chiral- 
odd transversity distribution I24L l25l . the time-reversal 
odd (T-odd) distribution [H EHIIlH El El EH and 
the Collins E3| T-odd fragmentation functions. 

In the partonic description of SIDIS, distribution and 
fragmentation functions depend on the scaling variables x 
and z, respectively (see below for the definition). At fixed 
and moderate values of the four momentum transfer Q 2 
and at large values of x and z, the contribution of multi- 
parton correlations or higher twist effects increases, even- 
tually leading to a breakdown of the partonic description. 
Model calculations indicate that SSAs are less sensitive 
to a number of higher order corrections than cross section 
measurements in both semi-inclusive 26] and in hard ex- 
clusive jH IH pion production. The measurement of 
spin asymmetries could therefore become a major tool 
for studyi ng q uark transverse momentum dependent dis- 
tributions 10, 

E1E1E1EHH and GPDslBlil in the 

Q 2 domain of a few GeV 2 . 

The beam-spin asymmetries in single-pion production 
are higher-twist by their nature [lj, |3(J, W$ and are ex- 
pected to increase at low Q 2 . Although large beam-spin 
asymmetries have been observed in measurements of ex- 
clusive electroproduction of photons (DVCS) [HHl], t ne 
only measurement of the beam-spin asymmetry in semi- 
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inclusive pion electroproduction was reported recently by 
the HERMES collaboration H at (z) w 0.4. Within sta- 
tistical uncertainties their value is consistent with zero. 

The cross section for single pion production by lon- 
gitudinally polarized leptons scattering from unpolarized 
protons may be written in terms of a set of response func- 
tions. The helicity (A e ) dependent part {(Jlu) 0, El 
arises from the anti-symmetric part of the hadronic ten- 
sor: 



dcr LU 



1 



The subscripts in olu specify the beam and target po- 
larizations, respectively [L stands for longitudinally po- 
larized and U for unpolarized). The azimuthal angle <j) is 
defined by a triple product: 



[fci x k 2 ] ■ P ± 
sm (b = — = — = — , 

\h x k 2 \\P±\ 

where k\ and k 2 are the initial and final electron mo- 
menta, and P_l is the transverse momentum of the ob- 
served hadron with respect to the virtual photon q. The 
structure function Ti' LT arises due to the interference 
of the longitudinal and transverse photon contributions. 
The kinematic variables x, y, and z are defined as: x = 
Q 2 /2(P iq ), y = (P^/iPih), z = (P 1 P)/(P 1 q), where 
Q 2 = —q 2 , q = k\ — k 2 is the 4-momentum of the virtual 
photon, Pi and P are the momenta of the target and the 
observed final-state hadron, and 7 2 = AM 2 x 2 y 2 /Q 2 . 

The beam-spin asymmetries in single-pion semi- 
inclusive leptoproduction were measured using a 4.3 GeV 
electron beam and the CEBAF Large Acceptance Spec- 
trometer (CLAS) [34| at JLab. Scattering of longitudi- 
nally polarized electrons off a liquid-hydrogen target was 
studied over a wide kinematic region. The beam polar- 
ization, frequently measured with a M0ller polarimeter, 
was on average 0.70 ± 0.03. Beam helicity was flipped 
every 30 msec to minimize the helicity correlated sys- 
tematics. The scattered electrons and pions were de- 
tected in CLAS. Electron candidates were selected by 
a hardware trigger using a coincidence between the gas 
Cherenkov counters and the lead-scintillator electromag- 
netic calorimeters. Pions in a momentum range of 1.2 
to 2.6 GeV were identified using momentum reconstruc- 
tion in the tracking system and the time-of-flight from 
the target to the timing scintillators. The total number 
of electron-7r + coincidences in the DIS range (Q 2 > 1 
GeV 2 , W 2 > 4 GeV 2 ) was «4x 10 5 . 

A critical issue in SIDIS processes is the assump- 
tion of factorization, i.e. that the hadron production 
cross section can be evaluated as a convolution of a in- 
dependent distribution function, a hard scattering and 
a z-dependent fragmentation function. This picture is 
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FIG. 1: Comparison of the distributions measured with 
CLAS at 4.3 GeV (circles) in x, Q 2 (GeV 2 ), missing mass 
Mx (GeV) , and the transverse pion momentum P± (GeV) with 
LUND-MC reconstructed events. The distributions are aver- 
ages over the range 0.5 < z < 0.8; the MC results are nor- 
malized to the same number of events. 



valid if the hadron originates from the fragmentation of 
the current quark, assuming there is sufficient energy so 
that the 4-momentum of the final hadron is not directly 
related to that of the struck quark. At low z hadrons 
may additionally originate from the fragmentation of the 
target, while at large z, in addition to higher twist ef- 
fects, diffractive effects are also important. Therefore a 
restricted range in 0.5 < z < 0.8 has been selected for 
the analysis. 

The event distributions in the restricted z-range have 
been compared with a Monte Carlo based on the LUND 
generator developed to describe high energy pro- 
cesses. In LUND, the pion production is dominated 
by direct production from string fragmentation, as op- 
posed to other processes such as target fragmentation and 
hadronic decays. Figure ^ shows comparisons between 
the experimental yields and the normalized MC distri- 
butions for different kinematical variables. The agree- 
ment of MC with the data suggests that the high energy- 
description of the SIDIS process can be extended to the 
moderate energies of this measurement. 

To verify the factorization ansatz, pion multiplicities 
have been extracted for different x ranges. Pion mul- 
tiplicities have been shown to be approximately equiv- 
alent to fragmentation functions |3j| which depend on 
the z variable only at fixed Q 2 . In Fig. [21 7r + multi- 
plicities normalized by the total number of events are 
shown as a function of z for different x bins. Within 
the range and the precision of the present measurement 
no x-dependence of multiplicities is seen. This experi- 
mental finding is also consistent with the assumption of 
factorization. 
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FIG. 2: Pion multiplicities as a function of z for different x 
ranges normalized by the total number of pions in the corre- 
sponding x-range. 



FIG. 3: The beam-spin azimuthal asymmetry as a function 
of azimuthal angle <j>, measured in the range 2=0.5-0.8. 



The azimuthal distribution of the beam-spin asymme- 
try for 7r + : 



A(d>) 



LU 



1 N+ - N~ 
PN+ + N- 



(2) 



is shown in Fig. [3J Here N ± is the number of events for 
positive/negative helicities of the electron and P is the 
beam polarization. The data show a clear sin <fi modu- 
lation from which a sin</> moment of 0.038 ± 0.005(stat) 
can be derived. 

The same quantity can be formed by extracting mo- 
ments of the cross section for the two helicity states 
weighted by the corresponding ^-dependent functions. In 
this case the sin</> moment is given by: 



.4 



sin <f) 
LU 



2 N 

i=l 



(3) 



The two methods are identical for a full acceptance 
detector, but in practice have different sensitivities to 
acceptance effects. In Fig. 0] the comparison of A s ^f 
derived with the two methods is presented as a function 
of the missing mass Mx evaluated in the ep — > e'n + X 
reaction. As can be seen the results agree well with each 
other over the full Mx range. The moments A s ^j can be 
also computed for each helicity state and for an average of 
zero helicity which corresponds to an unpolarized beam. 
These data shown in Fig. 01 provide an additional test of 
the absence of spurious azimuthal asymmetries. All these 
results indicate that within the statistical uncertainties, 
the acceptance corrections are small and under control. 

Contributions to the systematic uncertainties arise 
from spin-dependent moments of the cross section cou- 
pling to corresponding moments in the acceptance to pro- 
duce corrections to the measured sin <b moment 41 . The 
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FIG. 4: The beam-spin azimuthal asymmetry as a function 
of missing mass Mx, in 7*p — > 7r + X extracted in the range 
0.5 < z < 0.8. Triangles up and down are the results for posi- 
tive and negative helicities, respectively, and the filled circles 
are for their average. Open circles show the measured A 8 ^ 
extracted as a sin cj> moment of the spin asymmetry. Open 
squares show the measured A B ^ for the sample averaged 
over the beam polarization. Data are slightly shifted in Mx 
for clarity. 



contribution to uncertainties due to the CLAS accep- 
tance in all relevant kinematic variables {x,y,z,Pj_,<j)) 
is evaluated to be less than 0.004 in average and less 
than 0.007 in all bins. The systematic uncertainties in 
the measurement of the beam polarization contribute 
at an even lower level (0.002). Possible particle mis- 
identification over the accessible kinematic range changes 
the observed SSA by less than 0.001. To minimize radia- 
tive corrections, a cut on the energy of the virtual photon 
relative to the incoming electron (y < 0.85) was imposed. 
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The estimated radiative corrections do not exceed a few 
percent of the value of the SSA 37], and give a minor 
contribution to the systematic uncertainty. Other sys- 
tematic uncertainties are negligible. 

As can be seen in Fig. the missing mass on the rem- 
nant system is mostly occurring in the nucleon resonance 
region. Despite this, the beam SSA doesn't show a de- 
pendence on any specific final state. A sizable increase of 
the SSA only appears in exclusive n + production where 
the missing mass corresponds to the nucleon mass. For 
this reason, two different cuts on the missing mass Afj 
have been applied in the final analysis. A first cut at 
Mx > f -f GeV was chosen to exclude the contribution 
of the exclusive n + production on the nucleon. A higher 
cut at Mx > 1.4 GeV was also considered to reduce, in 
addition, the contribution of semi-exclusive n + produc- 
tion with a recoiling A°-resonance. For Mx > 1.4 GeV 
multi-hadron production is the dominant mechanism and 
the possible contribution of higher nucleon resonances 
in the recoiling system can be interpreted in terms of 
quark-hadron duality. This seems to be supported by 
the smooth and almost flat behavior of the data shown 
in Fig. 0] 

A further check is shown in Fig. where the z- 
dependence of the beam SSA is presented for increasing 
values of the missing mass cut. Due to the large correla- 
tion between the z and the Mx variables, an increasing 
Mx cut drastically reduces the number of events with 
large z, leaving almost unchanged the beam SSA. This 
indicates that, within the present statistical uncertain- 
ties, the fraction z of the virtual photon energy carried 
by the pion, rather than the missing mass Mx, is the 
relevant variable in the scattering process. 

Consistency with the factorization assumption, which 
has been already shown in Fig. [5] for ir + multiplicities, 
has also been investigated for the observable under study. 
In Fig. the beam SSA is presented as a function of z 
for different x-bins. Its general behavior suggested by 
the curve does not exhibit any significant ^-dependence, 
which is also consistent with factorization in the chosen 
kinematic range. 

The dependence of beam SSA on the tt + transverse 
momentum, P± , is shown in Fig. A linear dependence 
of SSA on the Pi (with kinematic zero at P± = GeV) is 
expected, at least for the moderate range of P± [llj, |l2| ■ 

The beam spin asymmetries averaged over the two spin 
states as a function of x and z are plotted in Fig. |H1 and 
listed in Table HJ Table [H] shows the relevant variables for 
the different x and z bins. The beam SSA is positive for a 
positive electron helicity over the entire measured range. 
The data do not show a significant x-dependence while 
the asymmetry is strongly increasing at high z, where 
according to LUND MC, the probability of the detected 
pion to carry the struck quark is maximal. 

The size and the behavior of the asymmetry are very 
similar for the two cases of missing mass cuts. For the 
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FIG. 5: The beam-spin azimuthal asymmetry as a function 
of z extracted for different cuts on the missing mass M x (in 
GeV), M x > 1.1 (circles), M x > 1.2 (squares), M x > 1.3 
(triangles up) and Mx > 1.4 (triangles down). 
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FIG. 6: Beam SSA as a function of z for different x ranges 
(same as in Fig. |5J for Mx > 1.1 GeV. The curve is a simple 
fit to all data to show the general behavior of the asymmetry. 



case of hi gher Mx cut, the data are compared with a 
prediction p3| based on the Sivers mechanism as the dy- 
namical origin of the observable. Within this framework, 
the asymmetry is given by the convolution of the T-odd 
parton distribution with the twist-3 fragmentation 
function E(x) j^H^. The latter function is responsi- 
ble for the strong z-dependence of the asymmetry. De- 
spite the fact that the formalism is much better suited 
for higher energy reactions, the agreement in size and 
behavior with the data is reasonable. 

The CLAS preliminary data on the beam SSA have 
been also interpreted in terms of the Collins mechanism 
US EJ El for a nrst determination of the twist-3 
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FIG. 7: Beam SSA as a function of P x for M x > 1.1 GeV 
(filled circles) and Mx > 1.4 GeV (open circles). 



distribution function e(x) |4*oL Elj . The magnitude of 
their e(x) is also consistent with predictions using the 
chiral quark soliton model 0, [3 C3 • 

At lower z HERMES reported results consistent with 
zero beam-spin asymmetry 4]. CLAS asymmetries, ob- 
tained at higher z and lower Q 2 , should be kinematically 
enhanced with respect to the higher e nerg y HERMES 
data Q, as pointed out in Refs. pffll Isij ". so that no 
evident contradiction can be inferred between the two 
measurements. 

The data at lower Mx demonstrate that for this ob- 
servable the transition from the semi-inclusive to the 
semi-exclusive and to the exclusive domains is smooth. 
In addition these data may provide a new field for test- 
ing, in the final state, the quark-hadron duality, which 
has been proved to work in the initial state for other 
observables. like the spin-independent [4fi| and the spin- 
dependent J47J structure functions, down to low values of 
Q 2 . 

In conclusion, we have presented the first measurement 
of a non-zero beam-spin asymmetry in single 7r + inclusive 
electroproduction in the DIS kinematic region. The aver- 
age asymmetry is 0.038 ±0.005 ±0.003 for a missing mass 
M x >1.1 GeV and 0.037 ± 0.007 ± 0.004 for M x >1A 
GeV. The asymmetry shows a strong enhancement at 
large values of z while no significant ^-dependence is 
present within the measured range. Detailed experimen- 
tal and Monte Carlo studies have been performed show- 
ing no large violation of the factorization assumption for 
the process and suggesting that the partonic description 
may be applied in the kinematical range of the measure- 
ment. New data are expected from experiments utilizing 
a 6 GeV polarized beam, allowing a more precise test of 
the factorization ansatz and the investigation of the Q 2 
dependence of the asymmetries. 



FIG. 8: The beam-spin azimuthal asymmetry as a function 
of x (a), in the range 0.5 < z < 0.8, and as function of z (b), 
in the range 0.15 < x < 0.4 for M x > 1.1 GeV (filled circles). 
The error bars show the statistical uncertainty, and the band 
represents the systematic uncertainties. The empty circles are 
results for Mx > 1.4 and have been slightly shifted to make 
them more visible. The empty square shows the HERMES 
result Q, which is an average over the range z=0.2-0.7 and 
£=0.02-0.4. The curve is a theoretical prediction [3l| . 
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TABLE I: SSA: x and ^-dependence for M x > 1.1 GeV 
(upper table) and Mx > 1.4 GeV (lower table). 



X 


/|Sinc£ l A _i_ A 
ST-LJJ 111 *-*stat =E t-*syst 


z 


-^-Llj^ i A s tat i /S. S yst 


0.18 
0.24 
0.31 
0.37 


0.041 ± 0.011 ± 0.004 
0.034 ± 0.008 ± 0.003 
0.053 ± 0.009 ± 0.004 
0.026 ± 0.012 ± 0.005 


0.54 
0.61 
0.69 
0.76 


0.017 ± 0.007 ± 0.002 
0.049 ± 0.009 ± 0.004 
0.062 ± 0.011 ± 0.004 
0.063 ± 0.014 ± 0.005 


0.18 
0.24 
0.31 
0.37 


0.043 ± 0.013 ± 0.005 
0.033 ± 0.011 ± 0.004 
0.046 ± 0.014 ± 0.005 
0.039 ± 0.023 ± 0.007 


0.54 
0.61 
0.69 
0.76 


0.027 ± 0.009 ± 0.003 
0.047 ± 0.013 ± 0.005 
0.076 ± 0.024 ± 0.005 
0.067 ± 0.080 ± 0.007 



TABLE II: Average values of Q 2 (GeV 2 ), W (GeV), y, P± 
(GeV) and z/x in each of x and z bins for Mx > 1-1 GeV 
(upper table) and Mx > 1.4 GeV (lower table). 



X 


(Q 2 ) 


(W) 


(y) 


(Pi-) 


<*> 




z 


(Q 2 ) 


(W) 




(x) 


0.18 


l.i 


2.5 


0.75 


0.46 


0.61 




0.54 


1.46 


2.3 


0.43 


0.27 


0.24 


1.3 


2.3 


0.67 


0.42 


0.61 




0.61 


1.44 


2.3 


0.43 


0.27 


0.31 


1.6 


2.2 


0.66 


0.41 


0.61 




0.69 


1.44 


2.3 


0.42 


0.27 


0.37 


2.0 


2.2 


0.67 


0.39 


0.61 




0.77 


1.43 


2.3 


0.36 


0.27 


0.18 


1.1 


2.5 


0.76 


0.43 


0.58 




0.54 


1.44 


2.3 


0.38 


0.26 


0.24 


1.4 


2.3 


0.70 


0.36 


0.57 




0.61 


1.41 


2.4 


0.34 


0.25 


0.31 


1.7 


2.3 


0.69 


0.32 


0.56 




0.69 


1.37 


2.4 


0.30 


0.23 


0.37 


2.0 


2.2 


0.70 


0.28 


0.56 




0.77 


1.26 


2.5 


0.24 


0.20 



